Extension of the electron dissipation region in collisionless Hall MHD reconnection 
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This paper presents Sweet-Parker type scaling arguments in the context of hyper-resistive Hall magnetohyrdo- 
dynamics (MHD). Numerical experiments suggest that both cusp-like and modestly more extended geometries 
are realizable. However, the length of the electron dissipation region, which is taken as a parameter by several 
recent studies, is found to depend explicitly on the level of hyper-resistivity. Furthermore, although hyper- 
resistivity can produce more extended electron dissipation regions, the length of the region remains smaller than 
one ion skin depth for the largest values of hyper-resistivity considered here-significantly shorter than current 
sheets seen in many recent kinetic studies. The length of the electron dissipation region is found to depend 
on electron inertia as well, scaling like {jn^/miY^^ . However, the thickness of the region appears to scale 
similarly, so that the aspect ratio is at most very weakly dependent on (me/mi). 



I. INTRODUCTION 

The problem of fast magnetic reconnection in collisionless 
or weakly collisional plasmas has become a subject of great 
interest in recent years, in large part due to its relevance to 
impulsive phenomena such as magnetospheric substorms, so- 
lar and stellar flares, and the sawtooth crash in tokamaks. A 
key development has been the realization that non-ideal terms 
in the generalized Ohm's law (including electron inertia, elec- 
tron pressure gradient, and Hall terms) provide a pathway to 
the onset of fast reconnection. This subject has a history span- 
ning nearly two decades, with important contributions from 
both the fusion and space plasma communities (cf. ^ 0] and 
references therein). It has been suggested that the Hall term 
in particular plays a key role in localizing the electron dif- 
fusion region. Within the framework of Hall MHD (or two- 
fluid) studies, there has been consensus on this point, though 
the issue of reconnection scaling has remained controversial. 
Some studies of reconnection in Hall MHD systems have re- 
ported a "universal" normalized reconnection rate of ~ 0.1, 
independent of dissipation mechanism B flp lla ] , or system 
sizeQUH. However, other studies [[13, IliOilMHQl [ill 
have found a broad range of dependencies on plasma param- 
eters such as the ion skin and electron skin depths, and on 
boundary conditions, appearing to refute the claim of univer- 
sality. 

In an interesting sequence of studies based on fully kinetic 
simulations employing particle-in-cell (PIC) methods iflTiflsl 
[T9I1 . questions have been raised regarding the conclusion, ob- 
tained from Hall MHD studies, that the electron diffusion re- 
gion is highly localized. In these PIC studies the electron dif- 
fusion region is not so localized. Rather it is elongated in 
the outflow direction to lengths on the order of lOrf; ( where 
di = c/bjpi is the ion skin depth) for electron to ion mass 
ratios of the order of 1/100, when Hall MHD simulations 
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produce localized (< di length) electron current sheets (see 
e.g. ll20ll2lll ). This raises the following interesting question: 
are fully kinetic simulations necessary to produce elongated 
electron current sheets? More specifically, is it possible to re- 
alize elongated electron current sheets within the framework 
of Hall or extended MHD models by means of a generahzed 
Ohm's law that includes additional closure terms and param- 
eterizes physical processes that may have kinetic origin? If 
the answer to this question is yes, it may facilitate greatly the 
capability of global multi-fluid models to represent the effects 
of kinetic physics at small scales. 

A recent study by Chacon et al. ll22ll presented scaling argu- 
ments and simulation results which indicated that in 2D elec- 
tron MHD (EMHD) systems with electron viscosity (hyper- 
resistivity), the electron dissipation region can take on a wide 
range of geometries; the governing equations permit the as- 
pect ratio (width/length) of the electron dissipation region to 
be of order unity (cusp like) or much smaller, i.e., leading to 
an extended electron current sheet. In the present work, we 
re-derive Chacon's EMHD scaling in perhaps a more trans- 
parent way, via Sweet-Parker type arguments about quantities 
at the upstream and downstream edges of the electron dissipa- 
tion region. Similar approaches have been employed recently 
by Malyshkin ll23ll . and also by Uzdensky ||24l] . These scaling 
results are then benchmarked numerically by two-fluid simu- 
lations of island coalescence both with and without electron 
inertia. We emphasize that neither this study, nor any other to 
date has produced a general method for analytically predict- 
ing the length of the dissipation region, which is perhaps the 
most important factor in determining whether the reconnec- 
tion in a given system is fast or slow. We show, in particular, 
that the length of the electron diffusion region cannot be sim- 
ply assumed to be a parameter fixed only by external bound- 
ary conditions, but depends on the mechanism that breaks the 
frozen-in condition (electron viscosity in this case). 

The organization of this paper is as follows. In section HIl 
we describe the simulation model and equilibrium profiles. 
In section |III1 we derive our scaling arguments. In section 
IV we describe the diagnostics used in analyzing simulation 
results. In section|V]we present and discuss the results of the 
simulations included in this study. The main conclusions are 
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FIG. 1: (color online) Initial out-of-plane current, Jz{x, y) (orange 
scale), and magnetic field lines (black contours). 



Our algorithm employs fourth-order accurate spatial finite 
differencing and the time stepping scheme is a second-order 
accurate trapezoidal leapfrog ll25ll26ll . We consider a square, 
2D slab with physical dimensions L x L = 12. Sd^ x 12. Sd^ 
(so that L = 12.8 in normalized units) and periodic bound- 
ary conditions imposed at a; = ±L/2, and y ~ ±L/2. The 
simulation grid is n^; x Uy = 512 x 512, yielding grid scales 
of Ax ~ Ay = 0.025. Note that compared to the simple 
form of Ohm's Law in resistive MHD, the generalized Ohm's 
Law in this two-fluid system (i.e. Eq. (|3]l above), contains 
four additional non-ideal terms: the Hall term (oc J x B), 
the electron pressure gradient term (Vpe), the electron iner- 
tia terms, (which are manifested by the terms in the def- 
inition of B'), and the hyper-resistive term (oc V^J). Note 
that of these non-ideal terms, only the hyper-resistive term is 
dissipative. Although hyper-resistivity (also known as elec- 
tron viscosity) is generally negligibly small in coUisional plas- 
mas, it can be significant in a coUisionless plasma. For a 
more thorough discussion of hyper-resistivity see, for exam- 
ple, Refs. SmillllMIIllilllllll. We consider the 
case of zero resistivity, and the model described here is thus 
intended to apply to coUisionless plasmas. 



summarized in section IVll 



II. SIMULATION MODEL 



A. Model Equations 



Our simulations are based on the two-fluid equations in- 
cluding hyper-resistivity. These equations in normalized form 
are lil: 



n (ftV, + V, • WO = J X B - Vp, 



dtB' 



-V X E', 



(1) 



(2) 



E' = -V,: X B + - (J X B' - Vpe) - t?hV2J, (3) 



dtn + Vj • V?i = -nV • 



(4) 



B' 



where p,, = nTe , Pi = nTi , p = Pi + Pe 
(1 - dlV^) B, and Ve = V, - J/n , J = V x B. 

For simplicity we assume an isothermal equation of state 
for both electrons and ions (qualitatively similar to the adi- 
abatic case), and thus take Te and Ti to be constant. The 
normalizations of Eqns. ([T])-(|4li are based on constant refer- 
ence values of the density uq and the reconnecting compo- 
nent of the magnetic field B^q, and are given by (normalized 



physical units): t 
where dj^e = c/upi^e- 
B ^ B/i?o, V,,, 

and J 3 /{riQeVA)- 



— cBqI {rriic), x x/di, 
.,JVa, where Va = ujcidt = 



and uj^i e 



B. Initial Equilibrium 

We begin with a two-dimensional (2D) equilibrium simi- 
lar to that used by Longcope and Strauss (1994), shown in 
Fig ([T]i. The normalized magnetic field and density profiles in 
our initial equilibrium are given by: 



'Bx{x,y) = Bq 
'By{x,y) = -Bq 
n{x,y) = 1 + 



7rw\ / 2t:x 
cos , — j sm 

cos j sin 

1-^2 



2(T, + Te) ■ 



(5) 
(6) 
(7) 



The boundary conditions are periodic in both directions. As 
required by these boundary conditions, B is periodic under 
y ^ y + L and x x + L . Note that the normalization 
parameter Bq is equal to 2n/L = 27r/12.8 = 0.4908. The 
density profile is chosen to satisfy the total pressure balance 
condition, which in normalized form is given by 



1 



-B' 



constant. 



(8) 



Unless otherwise stated, we take the (constant) total temper- 
ature to be Ttot = Ti + Te — 1.0, or in physical units 
AirnoTtot / BxQ = 1, so that the plasma (3 has a (minimum) 
value of 11.45 where B^ = 0.4908 and n = 1.3795. The 
density reaches a maximum value of = 1.5 in the center 
of the simulation domain {x = 0,y = 0), where B = 0. 
Initially the current is carried entirely by the ions. To prevent 
energy buildup at the grid scale, the simulations include fourth 
order dissipation in the density and momentum equations of 
the form fi4W^ where /i4 ~ 5.1 • 10~^. To avoid physically 
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artificial effects that can arise from exact reflection symme- 
tries of the initial condition, a small amount of random noise 
is added to the magnetic field and ion current at the levels 

\B,na.\ « 10-^ |J,„.„J « 10-5. 

This initial equilibrium configuration contains four mag- 
netic islands. Diagonally adjacent islands contain (out-of- 
plane) currents of like sign and are therefore mutually at- 
tracted. The islands with currents of opposite sign are of 
course mutually repelled. At rest these attractions and re- 
pulsions balance so that the equilibrium is metastable: ei- 
ther the upper right and lower left islands can coalesce at the 
center of the simulation domain, or the upper left and lower 
right islands can do so. Due to the symmetry of this sys- 
tem, whichever pair coalesces at the center, the opposite pair 
of islands will be repelled from the center coalescing at one 
of the corners of the simulation box. The symmetry of this 
metastable state is broken by giving the system an initial in- 
compressible in-plane velocity perturbation of the form: 



Y,{x,y) = Vo 



/27ry\ f2Trx\ 
sin(^— jx + sm(^— jy 



(9) 



Unless otherwise stated, the magnitude of this initial velocity 
is Vq = 0.1 in normalized units. 



III. SCALING ARGUMENTS 

In systems where the Hall term is important, the non-ideal 
region takes on a two-scale structure: an outer region where 
the ions are demagnetized but the electrons are still frozen in, 
and an inner region where the electrons are demagnetized as 
well. It should be emphasized that in the outer part of the non- 
ideal region, although the ions are demagnetized, there is no 
dissipation occurring there. Here we will refer to this outer, 
non-dissipative zone as the "Hall region." The inner, dissipa- 
tive region, where neither species is magnetized, and where 
(in our case) electron viscosity is the dominant term, will be 
referred to as the "electron diffusion region." Let Se, and Le 
denote the thickness and length of the electron diffusion re- 
gion, respectively. 

In the following scaling arguments we neglect the effects 
of electron inertia, which is valid when the smallest scale of 
interest is much larger than the electron inertial length, d^. 
We also suppose that the spatial scale of the electron dissi- 
pation region is small sufficiently small compared to the ion 
inertial length, di, so that ion motion may be neglected. In 
other words, we consider an EMHD model. In the regime 
where the above conditions are met, the dominant terms on 
the right hand side of the generalized Ohm's Law (Eq. (|3]l) are 
the Hall (proportional to J x B) term, and the hyper-resistive 
term (proportional to ryjj J). Let us consider the how these 
two terms conspire to influence the geometry of the electron 
region. Note that the arguments made here are very similar to 
those made by Uzdensky in a recently submitted paper ll24ll . 

From this point on we will work in a co-ordinate system 
rotated 45 degrees clockwise {x (l/\/2)(a; — y)) with re- 
spect to the system coordinates described in Section ( IIIBb . so 
that we call the inflow direction "y" and the outflow direction 



"x". In the case under consideration, there are two forces act- 
ing on the electrons as they flow away from the x-point: the 
J X B force accelerates the electrons away from the x-point 
along the x and z directions in a channel of width ~ 6e, while 
the hyper-resistive (electron viscosity) term resists the shear 
in the electron fluid caused by the locaUzed outflow. The peak 
electron outflow in steady state occurs at the point where these 
two opposing forces balance. This balancing of terms at the 
point of peak outflow is found to be well satisfied in our simu- 
lation results. Here we will take the half length (ie/2) of the 
layer to be the distance from the x-point/stagnation point to 
the location of peak electron outflow. Let B^y, Vex, and Jez 
denote the magnetic field, electron outflow velocity, and out- 
of-plane current density at the point {x,y) = {L^, 0). Then 
we have 



^ey J ez 

cue 

-^ey J ez 

cne 



0, 



neVp, 



VM- 



S'. 



u r x2 

9 9 : 



(10) 



where the approximation that ^ (5^^ holds provided that 
5e ^ Le- Between the inner and outer region, the electrons 
move at the E x B drift speed. So, at the downstream edge 
of the electron region, B^y = cEz/Vex- In steady state, 
i9tB = V X E = implies that Ez is spatially uniform. 

Thus, EziLe,o) = Ez{o,o) = my^Mo,o) ~ mBex/Sl 

(where Bex is the magnitude of the reconnecting component 
of the magnetic field at {x, y) = (0, ^e)) yielding: 



B. 



ey 



Vex5l ■ 



(11) 



Substituting this value for Bey into Eq. ( fTOl ) gives \2/ 



\Vex\ = \Vez\='^- -'^^'4^=^-4^, (12) 



ne ne 



inSe 



where Va = Bex/V^mTMn, Vae = Bex/V^Tnn^, d^ = 
c/ujpi, and de = c/uj^e- We must emphasize here the quan- 
tity {diVA) — {deVAe) — {cB ex / 4:TTne) is independent of 
the mass of either species. The present scaling arguments are 
essentially EMHD arguments, and as such there is no finite 
di. In the absence of electron inertia, there is also no de in 
an EMHD model. Uniformity of Ez holds across the electron 
current sheet as well, requiring 



'eyl 



Si 



(13) 



Finally, taking the plasma to be incompressible (rigorously 
valid in the EMHD limit), requires that \ Vey\Le — \Vex\Se, we 
can solve for Se, \ Vex\, \ Vey\, and Ez in terms of Va, di, ryjy, 
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and Le'. 



\Vey\ 
IK. I 

cE, 



{d^VA) 



1/3 



{d,VA)B,^ = 1_ (2Mk\ 



(14) 
(15) 
(16) 
(17) 



Equation ( fT4b predicts exactly the same scaling as that found 
by Chacon et al. (compare to Eq. (8) of reference |220). How- 
ever, as in that study, the length of the electron dissipation 
region remains as a free parameter This scaling permits (as 

1/2 

Stated by Chacon) both cases in which 5^ ^ ^ rj^^ , or 
5e ~ {riH /VaY^^ <C Le- However, this statement may lead 
one to believe that there exists a bifurcated solution, which 
permits only the two aforementioned scalings. In fact, neither 
the arguments presented here nor those made in Refs. |22] or 
ll23ll places any limit on the scaling of L^, thus failing to an- 
swer the question of whether the reconnection in any given 
system will be fast or slow. So in fact, a continuum of scal- 
ings may exist. The reconnection rate given in Eq. ( [Tt] ) is not 
explicitly dependent on r^n, but may implicitly depend on iih, 
depending on how L^. scales. Several limitations of the above 
scaling arguments are worth noting. First, the theory has as- 
sumed that ion motion is negligible on the spatial and tempo- 
ral scale of the reconnection process. Second, the effects of 
electron inertia have been neglected. Third, electron pressure 
gradients have been assumed to be negligible, largely for the 
sake of simplicity; on scales > di, electron pressure gradients 
are generally expected to be of the same order as the J x B 
term, as can be seen from Eq. ([3]). The first two of these points 
limit the applicability of these arguments to scales I that sat- 
isfy the condition de £ di. A further, more subtle point, 
which was mentioned above, is that the downstream force bal- 
ance criterion may not hold in all cases. However, even if 
there is no downstream force balance, the same scalings can 
be derived by enforcing steady state of the 2D EMHD equa- 
tion. In two dimensions, the magnetic field can be expressed 
without loss of generality as B = z x Wip + Bz{x, y)z. In this 
case, the in plane and out-of -plane components of the EMHD 
equation can be expressed (in normalized form) as: 



dti' = -B_L • V^B^ - ■r]H^±Jz = -E^ = const. 
dtB, = -Bj. • J, + mi^lfB, = 0. 



In steady state, again taking dx - 
(JjT^, these two equations become: 



(18) 




(19) 



and eliminating B^^ from these two equations, again yields: 
5e = {rinLe/ BxY/'^ , the same scaling given by Eq. ( fT4b . and 
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FIG. 2: (color online) Dissipation region in a simulation with 
rrie — 0. top panel: ID cuts of electron and ion outflow veloci- 
ties. Vex, & Vix- The vertical lines mark the location of peak elec- 
tron outflow (ends of electron diffusion region) and the points where 
electron and ion velocities come together, second panel: J^ix, y)\ 
dashed lines mark edges of electron diffusion region, third panel: 
Jiz {x,y); dashed lines mark the locations where the ion and electron 
flows come together, fourth panel: ion oulfiow,Vix{x,y), bottom 
panel: electron outflow. Vex {x, y) 



equivalent scalings come about for the other quantities as well. 
This approach is perhaps more general, as it doesn't require 
force balance at the downstream edge, though it may be less 
physically transparent. A general theory which predicts the 
length of the dissipation region in collisionless Hall MHD is 
in many ways the ultimate pursuit of all theoretical work on 
reconnection today, and as of yet no theory has delivered that 
result, except in special cases (e. g. dill). Lacking such 
a general theory, we now to numerical experiments to explore 
the dependence of on i]h- 
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rv. DESCRIPTION OF DIAGNOSTICS 



Before describing the results of our numerical experiments, 
we must describe the diagnostics employed in quantifying the 
extent of the non-ideal regions for each species. Near the cen- 
ter of the current sheet, the out-of-plane current density, J^, 
is due primarily to the electrons. In this study the width of 
the electron dissipation region is taken to be the full width 
at half max of the out-of-plane current, J^. The downstream 
edge of the electron dissipation region is taken to be the lo- 
cation of the peak in the electron outflow speed. Outside of 
the non-ideal region, the two species flow roughly together at 
the E X B drift velocity. Both the upstream and downstream 
edges of the non-ideal dissipation region are taken to be the 
first point outside the electron dissipation region at which the 
electron and ion velocities differ by less than, say, ten per- 
cent. For example, to determine the location of the down- 
stream edge, we begin at the downstream edge of the elec- 
tron dissipation region and search along x until the quantity 

1 — ^'l^^'gj ^ < 0.1. Figure ^ shows the shape of the ion 
and electron dissipation regions at late time in a simulation 
with = 0. The top panel shows ID cuts of Vex (solid), and 
Vix (dotted) along the outflow (x) direction. The inner pair 
of dashed vertical lines in this panel indicate the locations of 
peak electron outflow, which are taken to be the ends of the 
electron dissipation region. The outer pair of vertical lines 

mark the point where ( 1 — ^4^4^ I = 0.1, i.e. the ends of 



the Hall region. The second panel from the top shows the out 
of plane current, Jz{x, y), along with dashed lines indicating 
the extent of the electron dissipation region. The third panel 
shows Jiz, the out-of-plane ion current, with the dashed lines 
now indicating the extent of the Hall region in the x-y plane. 
The lower two panels show Vix, and Vex in the x-y plane. 
For comparison, corresponding results from a simulation with 
me/nii ~ 1/25 are shown in Fig. (O. The point of peak 
electron outflow is further downstream than that observed in 
the simulation with no electron mass. This extension of the 
current sheet was seen in an earlier study fxh, but was not 
well explored. Note that the electron outflow jets extend well 
beyond the boundaries of the electron diffusion region, and 
into the boundaries of the non-dissipative. Hall region. This 
feature is consistent with results PIC simulations. A recent 
paper by Hesse et al. ll36ll provides a detailed discussion of 
the non-dissipative nature of these outflow jets in the context 
of a kinetic system. The shape and extent of the Hall region 
are basically unchanged by varying r/H and nie/mi. So we 
will focus here on the electron dissipation region. In examin- 
ing the scaling geometry of the electron dissipation region, we 
will focus on quantities defined at the centers of the upstream 
and downstream edges. The upstream field, B^x, upstream 
density, Ue, and electron inflow speed Vey, are measured at 
the point {x, y) = (0, ^e), while the electron outflow speed is 
measured at the point (a;, y) = {Le, 0), which, as mentioned 
above is taken to be the location of the peak electron outflow. 
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FIG. 3: (color online) Dissipation region in a simulation with 
rrie/mi = 1/25. top panel: ID cuts of electron and ion outflow 
velocities, Vex, & Vix- The vertical lines mark the location of peak 
electron outflow (ends of electron diffusion region) and the points 
where electron and ion velocities come together, second panel: 
Jz (x, y) \ dashed lines mark edges of electron diffusion region, third 
panel: Jiz{x, y); dashed lines mark the locations where the ion and 
electron flows come together, fourth panel: ion outflow, Via; (a;, y), 
bottom panel: electron outflow, Vex {x, y) 



V. NUMERICAL RESULTS 

A. Hyperresistive Scaling 

Because, as mentioned above, no ab initio theoretical for- 
mulation of reconnection has yet successfully predicted the 
length of the cuiTent sheet in a Hall MHD system, we here 
begin our numerical investigation by describing the observed 
dependence of Le on varying hyper-resistivity, rj^- Figure (HI 
presents measurements of the length of the electron dissipa- 
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FIG. 4: (color online) Time series of the length of the electron diffu- 
sion region, Le, for varying levels of -qH (upper panel), and scaling 
of median plateau values of Le vs. rjH (lower panel). 



FIG. 5: (color online) Time series of the thickness of the electron dif- 
fusion region, Se, for varying levels of rjH (upper panel), and scaling 
of median plateau values of Se vs. rjH (lower panel). 



tion region from five two-fluid simulations of magnetic is- 
land coalescence. Time series of Lg are shown in the top 
panel of the figure. The five, solid curves represent the results 
of simulations with varying levels of hyper-resistivity, which 
do not include the effects of electron inertia. Notice that by 
tilci ^6.5 each of these simulations has achieved a relatively 
stable value of L^- The lower panel of the figure plots the 
median value of Le from the plateau phase (tild > 6.5) ver- 
sus rjH- The scaling of is quite clear, as shown by the 

1/3 

dashed line (of slope 77^ ) in the lower panel. Empirically, 

1/3 

Le oc 77^ in this modestly sized, periodic, coalescing sys- 
tem. How does this empirical result for Le relate to the scal- 
ings predicted in the previous section? If one simply inserts 
Le OC riK^ into Eq. ( fT4] i. one finds Se oc rj^^ . Time series 
of 5e are plotted in the upper panel of Fig. (|5]l. As was the 
case with Le, the thickness of the electron dissipation region, 
reaches a stable value after tild ^ 6.5. In the lower panel of 
Fig. the median plateau value of 6e is plotted versus rjH- 
The dotted line, plotted for reference, represents a scaling of 
77^^, in reasonable agreement with the data points. In the case 
of the lowest value of 77^ presented in this figure, the current 
sheet thickness drops nearly to the grid scale. So that point 
may not be entirely reliable. As a comparison to the results of 
Chacon et al., in Fig. (|6]l we plot 5e{Bex/ LeY^^, analogous 
to Chacon's Fig. 3. The upper dotted line represents 77]/'^, 

though the slope is actually closer to r/^^ (as shown by the 

1/3 

lower dotted line). Note that if Le ^ rjj^' (purely an empir- 
ical result), and 5e ^ 77^^ (in accordance with the predicted 
scaling of 5e, and the observed scaling of Le), then the aspect 
ratio of the dissipation region has an extremely weak depen- 
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FIG. 6: Scaling of 5e(Be./Le)'/^ vs. 77H 



dence on hyper-resistivity 5e/ Le ^ vli^ ■ Figure © presents 
the scaling of the aspect ratio of the electron dissipation re- 
gion. Time series of Se/Le are plotted in the upper panel of 
the figure. In the lower panel of the figure, diamonds repre- 
sent the median plateau value of the aspect ratio. As expected 
, the aspect ratio varies extremely little across the range of 
hyper resistivities considered here. These scalings of 6e and 
Le have been further tested in analogous systems of twice the 
physical size (L^ — Ly — 25.6) and half the physical size 
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FIG. 7: (color online)Time series of the aspect ratio of the electron 
diffusion region, 5e/Le, for varying levels of rjH (upper panel), and 
scaling of median plateau values of Se/Ls vs. rjH (lower panel). 

(Lj. = Ly = 6.4) of the system described above. The results 
are shown in Fig. (O. Not only the scaling of L^, but even 
its absolute numerical value is practically unchanged by dou- 
bling the system size, as shown by the triangles and diamonds 
in the upper panel of Fig. dS). The scaling of Le does how- 
ever break down if the system size is too small, in which case 
the small size of the coalescing flux bundles limits the poten- 
tial extent of Le- The dependence of Se on rjn, is unchanged 
by doubling or halving the system size. Thus for sufficiently 
large systems, the scalings found here do not appear to depend 
on the system size. This behavior appears to be quite differ- 
ent from that found in Ref. lfl6ll . which used a different initial 
condition and did not achieve steady state at any point during 
the evolution of the system. Next we turn our attention to the 
scaling of the reconnection rate. Time series of the electric 
field measurements are noisy. So, we consider measurements 
averaged over a window in time; 

where is the component of the magnetic vector poten- 
tial in the ignorable direction (i.e. the flux function up to a 
sign), and ( ) denotes an average over an interval of one ion 
cyclotron time. Here the reconnection has been normalized 
by the the local upstream values of the magnetic field and 
density, Bex, and rig, which as mentioned earlier are mea- 
sured at the point upstream where the out of plane current 
has fallen to half of its maximum value. Figure (|9|l shows 
normalized reconnection rates computed from simulation re- 
sults. Again the upper panel shows time series, while the 
lower panel shows median plateau values. These simulations 
achieve a stable plateau value for between t^lci — 6.5 



and trtci = 8.5. In the lower panel of Fig. (|9]) the diamonds 
represent the plateau value of normalized reconnection rate 
from the time series shown in the upper panel. The stars rep- 
resent plateau values of the reconnection rate over the same 
range of hyper-resistivities in the larger (Lx ~ Ly — 25.6di) 
system. The dotted line, plotted for reference represents a 
slope of rj^^^, as predicted by Eq. ([TtI i for our observed scal- 

1/3 

ing of Le oa r]fj . The observed reconnection rates appear to 
be within a factor of two of the predicted scaling. We note that 
although it is common practice to normalize the reconnection 
rate using local values of the magnetic field and density, these 
values are not necessarily values that fairly characterize the 
system. In this system for example, the meaningful physical 
purpose of the normalized reconnection rate is to character- 
ize the time scale of the process of reconnecting a significant 
amount of the flux that was contained in the islands at t = 0. 
To this end, it makes sense to use typical "global" values of 
magnetic field and density, which characterize the initial con- 
dition. Figure (fTOl i shows reconnection rates from the same 
simulations as in Fig. (|9]l. However, in Fig. ( fTOl i the reconnec- 
tion rates have been normalized to the global peak initial val- 
ues of B and n. Again the top panel shows time series, while 
the lower panel illustrates the scaling of the plateau value. 
These "globally" normalized reconnection rates appear to be 
nearly constant with a value of approximately 0.1 — 0.2 across 
the range of hyper-resistivities considered here. Note that the 
"raw," unnormalized reconnection rates are also nearly inde- 
pendent of while the upstream field and density do appear 
to depend on the dissipation coefficient. Figure (fTOl ) makes 
clear that although the locally normalized reconnection rate 
depends on rjH, the time needed to reconnected a significant 
fraction of the initial flux in this system (as reflected by the 
"globally" normalized reconnection rate) is roughly indepen- 
dent of the dissipation region, and in this sense, the reconnec- 
tion observed here is fast. 



B. Electron Mass Scaling 

We have observed above that the dimensions of the elec- 
tron current sheet can be extended in cases which include 
electron inertia, as compared to hyper-resistive simulations 
that do not include the effects of finite electron mass (see 
Figs.|4]and|5]). Additionally, for a given mass ratio, there ap- 
pears to be a threshold value of i]h below which Le is deter- 
mined by the mass ratio, and the effects of i]h become neg- 
ligible. The present study will not treat the general scaling 
of systems in which the effects of hyper-resistivity and elec- 
tron inertia are comparable. However, we briefly present scal- 
ing results from simulations with varying electron mass, for 
1] = 0,riH = 1.0 X 10~^. The top panel of Fig. [TT] shows 
time series of Le for mass ratios varying from me/m.i = 

and again the bottom panel plots the plateau values of 
Le VS. me/rrii. The dotted line in the lower panel represents 
an estimated scaling of Le oc {me/rrii)'^^^. This scaling is 
consistent with that found by Shay et al. lfl9ll , in the context 
of PIC simulations. The X in the figure denotes an extrapo- 
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FIG. 8: (color online)Scaling of the Length of the electron diffusion 
region, Le vs. hyperresistivity, r]H for varying system size(upper 
panel), and scaling of the electron diffusion region thickness, 5e 
(lower panel). 
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FIG. 9: (color online)Time series of the locally normalized reconnec- 
tion rate, {Er)l, for varying levels of rjn (upper panel), and scaling 
of median plateau values of (Er) l vs. rjH (lower panel) 
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lated value of Le « 15de (or about 1/3 of an ion skin depth) 
for a physical electron/ion mass ratio of 1/1836. Note that 
Shay et al. extrapolate to a value of Lp ~ O.Qdi for a phys- 
ical mass ratio, less than a factor of two different from our 
two-fluid result. This scaling can be understood using argu- 
ments presented in the aforementioned study by Shay et al.; 
the method is similar to that used earlier in this paper to de- 
rive hyper-resistive scalings. Suppose that the dominant terms 
in the electron momentum equation at the downstream edge 
of the electron dissipation region are the electron inertia term 
and the Hall (J x B) term. Then in steady state {dt = 0), and 
assuming the current to be entirely carried by the electrons 
at the scale of the dissipation region (the EMHD approxima- 
tion) we have the following steady state electron momentum 
equation: 



nmeVe • We — ne(V e X B), 



(21) 



Vh 



FIG. 10: (color online)Time series of the globally normalized re- 
connection rate, {Er)g, for varying levels of rjH (upper panel), and 
scaling of median plateau values of {Eaja vs. rjH (lower panel) 



whose X component can be expressed as: 



d_ 

dx 



r2 
' ex 



eVezBy. 



(22) 



This equation can be integrated from the x-point to the down- 
stream edge to arrive at an equation for Le- However, in order 
to integrate the above equation, we need to know something 
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FIG. 1 1 : (color online)Time series of the length of the electron dif- 
fusion region, Le, for varying values of me/rrii (upper panel), and 
scaling of median plateau values of Le vs. me/mi (lower panel). 



FIG. 12: (color online)Profiles of By{x) along the outflow direction 
at f — 11.0 in simulations with varying mass ratios. Note that the 
slope of By appears to be nearly independent of me/rrii at x=0. 
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about the x dependence of By. In agreement with the PIC re- 
sults of Shay et al. we have found that in our two-fluid 
model the profile of By along the outflow direction in the dis- 
sipation region does not appear to depend noticeably on the 
mass ratio, as can be seen in Fig. (fTZt . Taking Vez ^ Vac, 
assuming By to vary linearly with distance from the x-point 
as By{x) = B'yX, and integrating the above equation yields: 



3/8 



BoVa 



1/2 



Bo 

d^B'y 



3/4 



(23) 



This expression differs slightly from that obtained by Shay 
et al. They give the exponent of the (Bo/diBy) factor as 1, 
whereas we find it to be 3/ 4. Note that due to this difference 
our expression is entirely independent of rrii (as it should be 
since there was no nii in Eq. (l2Tl l. whereas theirs contains 
an implicit dependence (via the Va in the second factor). 
This difference is not likely to become apparent in numerical 
experiments, as neither Bq nor B'^ appear to depend on the 
mass ratio. Biskamp et al. ll39ll have made different scaling 
arguments using the observed uniformity of the canonical mo- 
mentum (F = jz — dlip) along the upstream edge of the diffu- 
sion region, which yield a scaling of Le ~ (me/mi)^/^. This 
scaling could also agree reasonably well with the results of 
our numerical experiment. Resolving the difference between 
a scaling of 3/8 and 1/3 is difficult with only one decade of 
mass ratios. However, it is encouraging that our two-fluid 
scalings reasonably reproduce the mass ratio scaling seen in 
PIC simulations. 

The scaling of the dissipation region thickness, Se, with 
me /mi is shown in Fig . [T3] The dotted line in the lower panel 
indicates a scaling of (me/m^)^/^. This is somewhat stronger 
than the scaling [Se/di ~ (me/m^)^/'*] observed by Daugton 
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FIG. 13: (color online)Time series of the thickness of the electron 
diffusion region, 5e, for varying values of me/rrii (upper panel), and 
scaling of median plateau values of 5e vs. me/rrn (lower panel). 



et al. lllTll . while showing good agreement with that observed 
by Dorfman et al. [Se/di ^ {me/piiV^'^] in kinetic simula- 
tions of reconnection in MRX pSl) (although in those simu- 
lations the frozen-in condition is broken by electron pressure 
agyrotropy rather than hyper-resistivity). Here the extrapo- 
lated valued of Se for a physical mass ratio is approximately 
Ide- The aspect ratio of the electron diffusion region does not 
appear to depend significantly on the mass ratio, since both 
the length and the width of the region are observed to scale 
with very similar powers of (nie/mi). 



VI. SUMMARY 

In this paper, we have made scaling arguments based on 
the balancing of terms in a hyper-resistive Ohm's law at the 
boundaries of the electron diffusion region, and have found 
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a scaling for 6p, equivalent to that found by other recent 
studies ll22l l23i |3 via a somewhat different method. How- 
ever, we emphasize here that these scaling arguments are ul- 
timately incomplete unless they make a prediction about the 
length of the dissipation region. The expression derived for 
the reconnection rate (Eq.[T7]i appears to be formally indepen- 
dent of electron mass, and hyperresistivity. However, formal 
independence is no guarantor of fast reconnection, because the 
upstream field and the length of the electron diffusion region 
may in fact depend on the dissipation mechanism. If is of 
order di or even smaller then reconnection rate can take on 
a numerically high value, but numerical results indicate very 
clearly that Le is a function of ijh for the simple coalescing 
system under consideration here. The reconnection rate nor- 
malized to local values of magnetic field and density depends 
on -qn in this system, which might lead one to conclude that 
the reconnection is dissipation dependent. However, when the 
reconnection rate is normalized to typical global values of B 
and n, it appears that the time scale of the reconnection pro- 
cess is in fact nearly independent of tjh, and could therefore 
be termed "fast". For the purposes of comparison with other 
systems, these global values in general yield a more mean- 
ingful value for the normalized reconnection rate, since local 
values may themselves depend on the dissipation mechanism. 

Electron inertia does affect the extent of the electron diffu- 
sion region. However, both the length and width of the elec- 
tron diffusion region appear to scale in very similar ways, so 
that the aspect ratio of the region appears to be independent 
of the mass ratio. The mass ratio dependence of the length 
of the electron dissipation region found here is in good agree- 
ment with PIC results. The present study has made no attempt 
to test the scaling of reconnection with varying guide field, 
Bz, or boundary conditions, which are left for future work. 
Additionally, we have considered electron inertia and hyper- 



resistivity separately in our numerical experiments. However, 
in the regime where both of these effects are significant, the 
parametric dependence of the reconnection region on both pa- 
rameters could be complicated. 

Recent PIC simulations of reconnection have produced 
highly elongated current sheets, which break up into many 
small islands. High Lundquist number simulations in the con- 
text of resistive MHD have also produced very extended cur- 
rent sheets, which are subject to a "secondary tearing insta- 
bility." This secondary instability also leads to the formation 
of many small plasmoids, similar to what is seen in the ki- 
netic systems. The present study was largely motivated by 
the question of whether these kinds of extended electron cur- 
rent sheets could be realized within the framework of Hall 
MHD. Electron viscosity does not on its own appear to be 
capable of producing highly extended current sheets. In col- 
lisionless kinetic systems, it has been shown that agyrotropic 
(off-diagonal) pressure tensor effects, which are not generally 
included in fluid models, play a dominant role in breaking 
the frozen-in condition. Hyper-resistivity models some of the 
physics associated with the electron pressure tensor. It is pos- 
sible that a more complicated fluid closure parameterizing the 
physics of the agyrotropic pressure tensor components may be 
capable of reproducing the current sheet extension and subse- 
quent plasmoid formation witnessed in kinetic systems. That 
investigation is a topic of presently ongoing research. 
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